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Abstract 

The radial variation of hot gas density and tem- 
perature in bright elliptical galaxies provided by x- 
ray observations can be used to accurately determine 
the radial distribution of total galactic mass using the 
condition of hydrostatic equilibrium. However, we 
set out here to solve the inverse problem. Starting 
with the known distributions of total mass and mass- 
losing stars in a specific large elliptical, NGC 4472, 
we attempt to solve the gas dynamical equations to 
recover the currently observed radial distribution of 
density and temperature in the hot interstellar gas. 
The galaxy is assumed to be initially gas-free as a 
result of early SNII-driven galactic winds. In seek- 
ing this agreement we consider a variety of assump- 
tions for the evolution of the hot interstellar gas: mass 
dropout, variation of stellar mass loss with galactic ra- 
dius, variation of the bolometric radiative cooling rate 
with galactic radius, and several supernova rates. Af- 
ter evolving for a Hubble time, none of these models 
accounts for several well-observed properties of the 
interstellar gas: gas temperatures that significantly 
exceed the mean stellar temperatures, positive tem- 
perature gradients within a few effective radii, shallow 
density gradients, and large interstellar gas masses. 

However, all of these discrepancies are lessened or 
disappear if large masses of hot gas are assumed to 
be present in the outer galactic potential at early 
times in galactic history. Most of the interstellar mass 
that contributes to cooling flows in galaxies like NGC 
4472 has not come from mass lost by galactic stars 
since ~ 1 Gyr. The sustained inflow of hot gas from 
this circumgalactic environment may help to resolve 
other long-standing problems that have beset models 



of galactic cooling flows: the wide range of L x /Lb 
for fixed Lb] the variable and often low iron abun- 
dances observed in ellipticals, and the failure (so far) 
to observe pronounced rotational flattening in x-ray 
images of slowly rotating giant ellipticals. 

The characteristic hot gas temperature profile ob- 
served in many bright ellipticals has a maximum at 
~3r c . This can be understood as the mixing of gas 
ejected from stars with old circumgalactic gas flowing 
in from the halo. Since the circumgalactic gas is hot- 
ter than the stars, mass dropout is not needed to flat- 
ten the gas density gradient. Moreover, the agreement 
of the stellar mass in NGC 4472 (and NGC 4649) 
within r e determined from stellar dynamics with the 
total mass determined by hot gas hydrostatic equilib- 
rium is an additional argument against mass dropout. 
Subject headings: galaxies: elliptical and lenticular - 
galaxies: cooling flows - X-rays: galaxies 

1. INTRODUCTION 

As the quality of x-ray observations of elliptical 
galaxies has increased, the modest agreement with 
simple theoretical models has diverged. In bright el- 
lipticals the x-ray emission is dominated by thermal 
radiation from optically thin interstellar gas. The ra- 
dial variation of density and temperature in the hot 
gas, p(r) and T(r), can be determined from observa- 
tions of the projected X-ray surface brightness and 
spectrum. But we show here that neither p(r) nor 
T(r) can be reproduced with gas dynamical calcula- 
tions based on variations of the standard assumptions 
of galactic cooling flow theory. 

In the usual theoretical approach the equations of 
gas dynamics are solved with source terms describing 
the rate that new gas flows into the interstellar gas, 
ejected by red giants in an old stellar population. The 
thermal energy of the interstellar gas is influenced by 
several effects: (i) the (kinetic) energy of newly in- 
jected gas that it inherits from the orbital motion of 
the parent stars, (ii) additional heating by Type la 
supernovae, and, most importantly, (iii) compressive 
heating in the galactic gravitational potential. The 
final gas temperature, T ~ 10 7 K, is comparable to 
the virial temperature of the galactic potential. Heat- 
ing by supernovae cannot be too large otherwise the 
iron abundance in the hot interstellar gas would ex- 
ceed observational limits or the heated gas would flow 
out of the galaxy as a wind, reducing the total x-ray 
luminosity L x below the threshold of detectability. 



Conversely, the interstellar gas cannot be free-falling 
rapidly in the galactic potential since the mass loss 
rate from evolving stars is insufficient to resupply the 
observed mass of interstellar gas in the frccfall time 
for large ellipticals, 10 9 yrs. For these reasons it 
follows that the hot interstellar gas is flowing very 
subsonically and the condition of hydrostatic equilib- 
rium is satisfied to an excellent approximation at most 
galactic radii. 

The main properties of gas in x-ray luminous el- 
lipticals have been confirmed with evolutionary gas 
dynamical calculations, but many important details 
remain unexplained. Some inconsistencies between 
theory and observations have been known for some 
time. For example: the wide scatter in x-ray lumi- 
nosity L x for galaxies of similar optical luminosity 
Lb (Eskridge et al. 1995) and the often surpris- 
ingly low (and varied) mean iron abundance in the 
hot gas compared to that observed in the stars (Ari- 
moto et al. 1997). It has also been recognized for 
many years that gas density slopes dp/dr (derived 
from surface brightness observations) are less steep 
than those predicted by the most straightforward the- 
oretical models. Recently, however, new inconsisten- 
cies have emerged. New observations of the gas tem- 
perature using ROSAT PSPC and ASCA (e.g. Davis 
and White 1996) indicate that the average gas tem- 
perature in massive ellipticals is ~ 1.5 — 2 times larger 
than the equivalent orbital temperature of the stars 
T* = {\im v jW)a1 where \i = 0.63 is the mean molec- 
ular weight and m v is the proton mass. The gas is 
much hotter than the stars. In addition gas tem- 
perature profiles T(r) in most or all of the massive 
ellipticals studied so far show an increase by 30 — 40 
percent from the central values to about three effec- 
tive radii where the temperature flattens or slowly de- 
creases toward larger radii (Jones et al. 1997; David 
et al. 1994; Trinchieri et al. 1994; Mushotzky et al. 
1994; Kim & Fabbiano 1995; Irwin & Sarazin 1996; 
Trinchieri et al. 1997). A plot showing these results 
can be found in Brighenti & Mathews (1997a). We 
show below that this peculiar temperature variation 
does not arise naturally in evolutionary flows based 
on gas lost from galactic stars. Finally, since all el- 
liptical galaxies rotate at some significant level, their 
x-ray contours are expected to be significantly elon- 
gated perpendicular to the axis of rotation out to an 
effective radius or more (Brighenti & Mathews 1996; 
1997b), yet most observed x-ray contours appear to 
be approximately circular! 



In view of the serious nature of these observa- 
tional inconsistencies with prevailing theory, in this 
paper we restrict our exploration of alternate theo- 
retical models to spherical geometry. Our approach 
is to begin with a simple gas-dynamical model for 
the evolution of hot gas in a realistic galactic poten- 
tial based on NGC 4472. We show that the final gas 
density and temperature variation in this model dis- 
agrees with observations of NGC 4472 in a variety of 
important ways. Then we vary the model parameters 
and assumptions to seek improvements, but instead 
only illustrate how insensitive the model is even to 
rather drastic alterations of the parameters. Finally, 
we show that the disagreement between the proper- 
ties of interstellar gas in our gas dynamical model and 
those observed in NGC 4472 can be alleviated if large 
amounts of additional gas are supplied near the be- 
ginning of the galactic evolution. The same desirable 
result may be achieved by adding gas to the galaxy 
more slowly over many Gyrs, but we do not explore 
this possibility here. The important point is that el- 
liptical galaxies can retain halos filled with very old 
hot gas. This gas may be a relic of an earlier time 
when these massive ellipticals were formed near the 
central focus of small groups. For some large ellipti- 
cals in Virgo this ancient gas has survived after the 
galaxy has entered the environment of this rich clus- 
ter. 

2. NGC 4472 AND THE STANDARD ISM 
MODEL 

Recently we have reviewed the x-ray observations 
of three bright nearby ellipticals: NGC 4472, 4636, 
and 4649. From the known gas density and tempera- 
ture profiles in these galaxies and their optical prop- 
erties, we determined the variation of total and stellar 
mass densities as functions of galactic radius assuming 
that the hot interstellar gas is in hydrostatic equilib- 
rium (Brighenti & Mathews 1997a). For two of these 
galaxies, NGC 4472 and 4649, the hot interstellar gas 
is in equilibrium with the potential of the stars for 
r < r e , the effective radius, but the potential of the 
massive dark halos dominates beyond r e . For these 
two galaxies the agreement of the total mass and den- 
sity with that of the stars alone over 0.1 < r/r e < 1 
verifies that the mass to light ratio determined from 
the stars and the gas temperature determined from 
x-ray spectra are accurate and that hydrostatic equi- 
librium maintains. We can therefore be confident that 
the underlying mass structure of these galaxies is well 
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understood. 

2.1. Properties of NGC 4472 

For a realistic galactic model in which to test 
various gas dynamical models described below, we 
use the total mass M tot (r) = Md a rk{r) + M*(r) of 
NGC 4472 found by assuming hydrostatic equilib- 
rium using gas temperatures and densities based on 
x-ray observations. The variation of gas density n(r) 
for NGC 4472 plotted in Figure la combines results 
from ROSAT HRI and PSPC observations of Irwin 
& Sarazin (1996) and Einstein HRI observations of 
Trinchicri, Fabbiano & Canizares (1986). The an- 
alytic fit to the gas density n(r) = p(r)/1.25m p in 
Figure la is sum of functions n(r) = £rii(r) where 
m(r) = n (i)[l + [r/r (i)] pW ] _1 . For 4472 we used 
three such functions with the following parameters: 
n (i) = 0.095, 0.00597, -.0004 cur 3 ; r (i) = 0.107, 
0.95, 10. in units of r e ; p(i) = 2.0, 1.14, 1.19. The 
gas temperature profile for NGC 4472 in Figure la is 
fit with a function of the form T(r) = 2T m [r m / (r + 
fot) + (f/rrn) 5 ] -1 w i tn parameters T m = 0.75 x 10 7 
K, r m = 0.5r e , r ot = 0.75r e and q = 0. Although 
the observed temperature is an average along the line 
of sight at projected radius R weighted by the local 
emissivity ~ p 2 , we assume that the temperature at 
physical radius r is the same, i.e. T(r) « T(R). This 
approximation is appropriate because of the steep de- 
cline of p with galactic radius (Brighenti & Mathews 
1997a). 

The total mass that confines the hot interstellar 
gas is found from the condition for hydrostatic equi- 
librium 
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M tot (r) is plotted in Figure lb. Also shown in Fig- 
ure lb is the de Vaucouleurs stellar mass distribu- 
tion M* (r) determined from the optical data for NGC 
4472: luminosity L B = 7.89 x 10 10 L B ,o, distance 
D = 17 Mpc, effective radius r e = 8.57 kpc, stel- 
lar mass to light ratio 9.20 (van der Marel 1991), and 
total stellar mass M» t = 7.26 x 10 11 M Q . M tot and 
the local total density (Figure lc) are dominated by 
the dark halo for r > r e but within r e the x-ray dis- 
tribution is determined by the stellar potential alone. 
As noted earlier, the total mass in 0.1r e ^ r Js r e is 
in excellent agreement with the expected stellar mass, 
M to t ~ M*. The galactic potential GM tot /r that we 



adopt for our dynamic models for gas flow in NGC 
4472 is based on M toi (r) shown in Figure lb. 

Iron abundance determinations for the hot gas in 
NGC 4472 vary with method and instrument used: 
z Fe = 0.20 (Scrlcmitsos et al. 1993 with BBXRT), 
z Fe = 1 - 2 (Forman ct al. 1993 with ROSAT), z Fe = 
0.63 (Awaki et al. 1994 with ASCA), and z Fe = 1.18 
(Buote & Fabian 1997 with ASCA). Some, but not all, 
of this variation is due to significant differences in the 
adopted abundance for the sun or solar system (Ishi- 
maru & Arimoto 1997). Arimoto et al. (1997) adopt a 
mean (emission- weighted) gas abundance z Fe — 0.33 
in solar units with z Fe Q = 4.68 x 10~ 5 . In our models 
much of the interstellar iron comes from stellar mass 
loss. Gas ejected from stars at any radius is assumed 
to have the same iron abundance as the local parent 
stars. For the stellar iron abundance in NGC 4472 we 
assume a solar value at the galactic center which de- 
creases as z %Fe oc r~ - 3 far from the galactic core (Ari- 
moto ct al. 1997), i.e. z* Fe (r) = 1.0/[1 + (r/r b ) 2 } - 15 
where r& = 200 pc is the core or "break" radius of 
NGC 4472 (Faber et al. 1997). We have deliberately 
chosen a somewhat smaller central stellar iron abun- 
dance than the z* Fe (0) = 1.6 adopted by Arimoto et 
al. (1997); previously it had been thought that z* Fe 
followed the spectral features of Mg in the solar ratio, 
but it is now thought that Mg/Fe exceeds the solar 
ratio in bright ellipticals (Scott Trager, private com- 
munication). Since all of these parameters are highly 
uncertain, the absolute iron abundances we calculate 
are very approximate. Our main interest here are the 
relative emissivity- weighted abundances {z Fe /z FeQ ) 
among the various dynamical models. 

2.2. The Standard Model for the Hot Gas in 
NGC 4472: Model STD 

The standard gas dynamics equations that describe 
the evolution of hot interstellar gas in ellipticals are: 
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Here e = 3kT/2pm p is the specific thermal energy. 
Both stars and supernova contribute mass to the in- 
terstellar gas, a = a* + a sn , but the contribution 
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from supcrnovae a sn is very small. An old stellar 
population expels gas at a rate a*(f)p* gm cm~ 3 s _1 
where a*(t) = a*(tn)(t/*ra) _1 ' 3 > *n = 13 Gyrs is the 
current time and a*(t„) = 5.4 x 10~ 20 s^" 1 (Math- 
ews 1989). The time variation of a*(t) is surprisingly 
insensitive to the IMF assumed, but a*{t n ) is only 
know to within a factor of ~ 2. As described above, 
M to t (r) is known directly from the observed gas tem- 
perature and density in NGC 4472 and the stellar 
density /?* (r) is based on a de Vaucouleurs profile nor- 
malized to values appropriate to NGC 4472. The loss 
of thermal energy by optically thin thermal emission 
is described by the cooling coefficient A(T) modeled 
after the results of Raymond, Cox & Smith (1976). 
The gas is assumed to be heated only by supernova 
of Type la since we assume that the early phase of 
galactic winds generated by Type II supernovae has 
already subsided at time t = 1 Gyr when we begin 
our calculations. The mean gas injection energy is 
e Q = 3kT /2/j,m p where T a = (a*T*+a an T sn )/a. The 
local stellar temperature T* (r) is found by solving the 
equation of stellar hydrodynamics (with isotropic ve- 
locity ellipsoids) in the total mass potential of NGC 
4472. Supernova heating of the interstellar gas is de- 
scribed by 

a sn T sn = 2.13 x 1(T 8 SNu(i) (£ sn /10 51 ergs) 

h- 1 - 7 (L B /L BQ )-°- 35 K s- 1 

where h = ff/100 = 0.75 is the reduced Hubble con- 
stant and we adopt E sn = 10 51 ergs as the typical 
energy released in a supernova event. This equa- 
tion incorporates the mean stellar mass to light ra- 
tio for elliptical galaxies from van der Marel (1991), 
M* t /L B = 2.98 x IO-^Lb/Lbo) - 35 /* 1 - 7 . (The in- 
dividually determined mass to light ratio for NGC 
4472, M* t /L B = 9.20, is slightly lower than that 
of the average elliptical having the same L B .) For 
our galaxy model (based on NGC 4472) we take 
L B = 7.89 x 1O 1O L B 0- Although there is some in- 
formation available about the current Type la super- 
nova rate in ellipticals (Turatto, Cappellaro & Benetti 
1994), almost nothing is known about supernova rates 
in the distant past. For this paper we adopt a power 
law to describe a varying rate of Type la supernovae: 

SNu(t) = SNu(t n )(t/t n )- p 

with p = 1 and SNu(t n ) = 0.04 SNu (1 SNu = 1 
supernova every 100 years in stars having luminosity 
L B = 10 10 Lbq). We assume each supernova ejects 



1.4 M© into the local interstellar gas and that half of 
this, 0.7 Mq, is iron. The iron abundance currently 
observed in the hot gas provides a global constraint 
on the past supernova rate and our parameters are 
chosen with this constraint in mind. Finally, the de 
Vaucouleurs model for the stellar component in our 
model galaxy extends out to radius r t = 100 kpc 
where the stellar density is very low, but the dark 
halo continues out to 470 kpc where its total mass is 
Mdark,t = 2.7 x 10 13 Mq. Since the gas is ejected from 
stars within 100 kpc, the much more extensive dark 
halo - which may resemble those of central ellipticals 
in groups - is irrelevant to most of the following gas 
dynamical solutions. 

The gas-dynamical equations are solved using a 
standard one-dimensional second order Eulerian code. 
We use reflection boundary conditions at the origin 
and outflow boundary conditions at the outer bound- 
ary, r = 470 kpc. 

2.3. Comparison of Standard Model (STD) 
and Observations 

Along with the observations of NGC 4472 in Figure 
2 we also show the current variation of density and 
temperature predicted by solving the gas dynamical 
equations above. These equations are solved begin- 
ning with no gas in the galaxy at t = 1 Gyr and car- 
ried forward to the current time, t n = 13 Gyr. This 
gas-free initial condition can be understood if SNII ex- 
plosions had driven strong winds out of the galaxy at 
early times (Mathews 1989; avid et al. 1991; Loewen- 
stcin & Mushotzky 1996). The elliptical in our gas 
dynamical model is assumed to be isolated so that 
gas can flow out beyond the galaxy if it wishes. 

The deficiencies of the model are apparent! The 
predicted density gradient for this "standard model" 
(STD) shown in Figure 2 is steeper than the observed 
density distribution throughout the entire galaxy. 
Traditionally this discrepancy has been recognized 
only in the central regions. The projected x-ray sur- 
face brightness of the model galaxy E x (i?) is also 
steeper than that observed for NGC 4472. This den- 
sity discrepancy still persists for models, not shown in 
Figure 2, in which the calculation is begun at t = 0.1 
Gyr to increase (by ~ 3) the total amount of gas 
ejected from the stars. Although the gas tempera- 
ture predicted by our model is roughly comparable 
with observed values shown in Figure 2, it is defi- 
nitely too low at all radii by a factor of ~ 2 — 3 
except perhaps within a few kpc from the center. 



4 



In addition, the shape of our computed temperature 
variation T(r) does not have the characteristic pos- 
itive gradient dT/dr out to r ~ 3r e that is typi- 
cal of most or all galactic cooling flows observed so 
far (Brighenti & Mathews 1997a). For deep gravita- 
tional potentials, as for example where the de Vau- 
couleurs mass dominates, the temperature gradients 
determined by our gas dynamical model are slightly 
negative at most radii. This negative gradient is due 
to the work done by gravitational compression; as gas 
expelled from stars at radius r loses energy by radi- 
ation, it sinks deeper in the galactic potential. For 
shallower potentials, however, dT/dr can be positive, 
but such potentials are not realistic models for galax- 
ies like NGC 4472. We have not illustrated the ve- 
locity field corresponding to the model in Figure 2 
since it is not directly observed. The flow velocity is 
negative throughout most of the galaxy but becomes 
slightly positive in the very tenuous gas at r > 160 
kpc. Within this radius the velocity gradually in- 
creases inward reaching -20, -50, and -140 km s _1 at 
10, 1, and 0.32 kpc respectively. 

The mean interstellar iron abundance in model 
STD, weighted by p 2 A, is {z Fe ) = 1.03 in solar units 
which is higher than the mean (mass-weighted) stel- 
lar value z- tFe « 0.34. If the observed interstellar 
abundance in NGC 4472 is as low as Arimoto et al. 
(1997) claim, z Fe m 0.33, then our (rather low) cur- 
rent supernova rate SNu(i„) = 0.04 may still be an 
overestimate. 

In model STD essentially all of the gas lost from 
the stars is retained within the galaxy. Most of the 
gas created in model STD from t = 1 to 13 Gyr, about 
7.2 x 10 10 Mq, eventually cools at the very center of 
the flow leaving only a relatively small mass of hot 
gas, 0.6 x 10 10 Mq, to produce the x-ray luminosity 
L x . In more realistic rotating cooling flows the gas 
cools to a disk (Brighenti & Mathews 1996; 1997b). 
In Table 1 we list several global parameters of the 
standard model STD and its variants described below. 

3. VARIATIONS ON THE STANDARD MODEL 
3.1. Mass Dropout: Model DO 

The conventional theoretical means of flattening 
steep density profiles in galactic cooling flow models 
is straightforward: remove hot gas from the cooling 
flow near the galactic center! In addition to correct- 
ing the density profile, this much-used assumption 
of "mass dropout" was also motivated by the pos- 



sibility that thermal instabilities could develop from 
hypothetical inhomogencities in the cooling flow gas 
(Fabian & Nulsen 1977; Stewart et al. 1984; White 
& Sarazin 1987a, 1987b, 1988; Thomas et al. 1987; 
Vcddcr, Trcster & Canizares 1988; Sarazin & Ashe 
1989; Bertin & Toniazzo 1995). However, such ther- 
mal instabilities do not grow in the expected man- 
ner. Low amplitude thermal instabilities are gen- 
erally stable (Balbus 1991) and coherent non-linear 
perturbations of larger amplitude oscillate radially in 
the cooling flow atmosphere (Loewenstein 1989) so 
that they are alternatively overdense and underdense 
relative to ambient gas, reducing or eliminating the 
growth rate of the instability. Even if instability were 
possible, such initially overdense, oscillating regions 
are not likely to remain coherent since hydrodynami- 
cal models show conclusively that they disintegrate 
by Rayleigh- Taylor and Kelvin-Helmholtz instabili- 
ties before completing a single oscillation (Hattori & 
Habe 1990; Malagoli et al. 1990). Moreover, the 
origin of the required inhomogeneities is also unclear 
since there is no known source of entropy perturba- 
tions on the appropriate scales in the interstellar gas 
(Mathews 1990). 

Nevertheless, the problem of steep theoretical den- 
sity profiles is clearly still present in Figure 2! In 
view of this and in deference to those many as- 
tronomers who have invoked mass dropout, we now 
describe some models for NGC 4472 using the stan- 
dard dropout assumptions which are expected to flat- 
ten dp I dr. When dropout is included, an additional 
term —qp/tdo is introduced on the right hand side of 
the equation of continuity but the equations of motion 
and thermal energy remain unchanged. The charac- 
teristic dropout time is usually assumed to be the 
time for local radiative cooling, tdo = 5m p kT/2p 1 pA 
(e. g. Sarazin & Ashe 1989), and the constant di- 
mensionless parameter q controls the amount of mass 
dropout. Gas is removed from the flow on the spot, 
i.e. tdo is assumed to be much less than the local flow 
time. Although tdo <C tfi ow w r/\u\ is widely used 
in mass dropout calculations, it is often strongly vi- 
olated. The rate that gas is removed from the flow 
is proportional to pt^ cx p 2 and therefore increases 
closer to the galactic center as required to flatten 
dp/dr. The temperature of the remaining gas that 
docs not participate in the dropout is increased. This 
occurs since after dropout less gas is present at small 
radii to support the weight of the same large global 
mass of interstellar gas so compressional heating is in- 
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tcnsificd at small radii to maintain hydrostatic pres- 
sure equilibrium. However, as gas cools and under- 
goes local dropout it will emit radiation from system- 
atically cooler gas. The net temperature observed at 
any radius is therefore a weighted mean between the 
background flow and the gas that is dropping out. 
The local x-ray emissivity cab into the ROSAT band 
(AE = 0.2-2.4 kcV) is increased by a factor (l+qA ) 
and the local mean temperature is lowered from T to 
T eff = T(l + qAi)/(l + qA ) where T is the temper- 
ature of the gas that has not yet dropped out and A 
and Ai are defined by 

A =— — [ T (TT—dT' 

A " " T*+" A AE J, {1) A 

The functions A (T) and Ai(T) are plotted in Figure 
3. 

Figure 4 shows the density and temperature struc- 
ture of a cooling flow with mass dropout using q = 1.2. 
A comparison with Figure 2 reveals that dp/dr is in- 
deed flattened within about 40 kpc by the dropout 
but the flow at larger radii, where the dropout is 
very small, is unchanged. The slope dp/dr between 
1 and 10 kpc is almost flattened to the slope of the 
NGC 4472 density profile and a further increase in q 
would bring the slope there into agreement. However, 
dropout has also lowered the gas density in the bright 
inner galaxy so a«(t„) would need to be increased by 
> 4 to bring the dropout solution into approximate 
agreement with local NGC 4472 densities within 10 
kpc. Such a large increase in a*(t n ) may be incon- 
sistent with reasonable stellar IMF slopes and mass 
cutoffs (Mathews 1989). In any case, the gas density 
is very much lower than observed values at r ~ 100 
kpc where the x-ray observations are still very reli- 
able. Although the beneficial influence of dropout on 
the density profile near the core has long been recog- 
nized, its failure to flatten dp/dr at larger radii has 
been unappreciated probably because most previous 
dropout solutions also used the steady state approxi- 
mation which is known to diverge and become unre- 
liable at large galactic radii (e.g. Vedder, Trester & 
Canizares 1988). 

But the most serious problem for the cooling flow 
with dropout is seen in the temperature profile in Fig- 
ure 4. The light solid line shows the temperature 
variation of the background gas which is increased 
over that in Figure 2 as explained above. When 
the cooler temperatures of locally dropping out gas 
is included, however, the apparent gas temperature 



T e f f (r) (heaver solid line in Figure 4) is slightly lower 
than the temperature without dropout shown in Fig- 
ure 2. According to the usual dropout assumptions, 
the dropout causes a lowering of the gas temperature 
even at very large galactic radii (Figure 4) since the 
small mass that drops out there is compensated by 
the much longer dropout time tdo- Therefore, even if 
mass dropout were a physically acceptable process, it 
fails to correct the deficiency of gas at large galactic 
radii (large negative dp/dr) and moves the computed 
gas temperatures further from observed values. 

Aside from the many theoretical arguments against 
the widely-used assumption of mass dropout, x-ray 
observations also indicate that galactic cooling flows 
cannot have a pronounced global multiphase char- 
acter in which regions of low entropy are cooling 
out. Unless there is a conspiracy of compensating 
effects, the near perfect agreement of the total and 
stellar mass in Figure lb - and also for NGC 4649 
(Brighcnti & Mathews 1997a) - are inconsistent with 
mass dropout. In the presence of dropout the ap- 
parent temperature is lowered from T to T e ff which 
reduces M to t ( r ) determined from the equation of hy- 
drostatic equilibrium. Also the density gradient is 
flatter, as in Figure 4, implying a smaller confining to- 
tal mass than actually exists. But the density shown 
in Figure 4 represents that of the background flow 
and has not been corrected for emission from regions 
of higher density that are dropping out; so the reduc- 
tion in M tot (r) is not as large as would appear from 
Figure 4. If mass dropout is really present and not 
allowed for in determining M tot (r) from the equation 
of hydrostatic equilibrium, Gunn & Thomas (1996) 
have shown that the apparent total confining mass 
is lowered by 20 to 60 percent. In typical dropout 
models, gas is removed from the flow instantaneously. 
But the pressure gradient in the flow readjusts only 
after a local sound crossing time, too slow to follow 
the effects of dropout on the gas density and tempera- 
ture. However, the total mass (Figure lb) determined 
from observed gas p(r) and T{r) (Figure la) agrees 
to within 10 percent of the independently determined 
stellar mass for 0.1r e < r < r e . This is strong ob- 
servational evidence that distributed mass dropout is 
not dominating the cooling flows observed in these 
two well-observed ellipticals. 

3.2. Non-uniform a*: Model NUA 

Since mass dropout has failed to account for the 
observed interstellar gas in NGC 4472, we abandon 
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caution and explore models in which the fundamen- 
tal coefficients in the source terms of the gas dy- 
namical equations are varied in an arbitrary fash- 
ion. For example, suppose proportionally more gas 
is ejected from stars further from the galactic core, 
then it might be expected that the ISM density pro- 
file would flatten. In an attempt to investigate this 
possibility, we impose an arbitrary spatial depen- 
dence on the stellar mass loss coefficient in which the 
stellar mass loss is enhanced at larger radii: a* = 
<x*{t n ) (r/r a ) a (i/t„)~ 1,3 . We have considered vari- 
ous values of a > and r a . A representative example 
of the resulting flow at t n with a = 0.4 and r a = r e is 
shown in Figure 5. Even though more (less) gas is be- 
ing expelled from the stars at large (small) radii, the 
final solution for both n(r) and T(r) is almost identi- 
cal to that in Figure 2 with uniform a.J The overall 
vertical normalization has changed slightly since the 
total mass lost from the stars M* t J a^dt is somewhat 
different, but the slope of the final density profile is 
very insensitive to the radial dependence of the stellar 
mass loss rate for all of the a-r a pairs we have consid- 
ered. Evidently this result is due to the insensitivity 
of the gas temperature profile to a* and the neces- 
sity for the gas to adjust to just that radial density 
distribution that achieves hydrostatic equilibrium. 

3.3. Non-uniform A: Model NUL 

Since none of the previous modified solutions has 
increased the gas temperature toward observed val- 
ues, particularly further out in the galaxy, perhaps 
an adjustment of the radiative cooling rate can ac- 
complish this. Although the total radiative cool- 
ing coefficient A(T) changes slowly with tempera- 
ture, A(T = 10 6 ) « 3A(T = 10 7 ), an additional 
dependence on metallicity should also be present at 
some level and this has not been considered in the 
models previously discussed. Since the metallicity 
z of the hot cooling flow gas increases inwards and 
dA/dz > 0, we consider the evolution of cooling 
flows with arbitrarily enhanced cooling at small radii: 
A(T, r) = (r/r s ) s A(T) where s < and A(T) is the 
cooling rate used in the standard solution (Figure 2). 
We have solved the cooling flow hydrodynamics with 
several s-r s pairs (keeping all other parameters iden- 
tical to those of the standard model) and the results 
for s = —0.6 and r s = 3r e are shown in Figure 6 at 
time t n = 13 Gyr. 

As a result of this rather large adjustment in A the 
central gas temperature has increased slightly at r < 



3r e « 26 kpc, but the gas temperature is still much 
lower than that observed in NGC 4472. However, at 
r > 3r e the computed temperature T(r) is almost 
unchanged from that of the standard model in Figure 
2. The slope of the density profile dn/dr in Figure 6 is 
significantly improved within about 10 kpc, but is still 
much too small further out. Our results for other s- 
r s pairs are similar so we must conclude that drastic, 
unphysical adjustments to the radiative cooling rate 
cannot bring the model temperatures (or densities) 
into satisfactory agreement with observations. 

3.4. Higher Supernova Rate: Models SN1 
and SN2 

We now explore the possibility of raising the gas 
temperature by increasing the supernova rate above 
that used in the standard solution of Figure 2. If the 
gas can be heated in this manner its pressure and den- 
sity scale heights will also increase, flattening dn/dr 
and improving solutions for both T(r) and n(r). Of 
course as the supernova rate SNu(t) is increased, we 
can anticipate that the iron abundance in the hot gas 
will also increase even further beyond the observed 
mean value in the hot gas (Loewenstein & Mathews 
1991). But we shall avert our concern about the iron 
abundance for the time being, hoping that a sepa- 
rate solution for reconciling the iron abundance can 
be found later by some exotic theoretical artifice. 

The supernova rate SNu(t) = SNu(i„)(t/t n )" p 
(where p > 0) can be altered either by raising p to 
increase the past rate of Type la supernovae or by in- 
creasing the coefficient SNu(t„) which increases SNu 
at all times. Our value p = 1 for the standard solu- 
tion allows for some reasonable cosmic evolution but 
is less than the exponent 1.3 in the stellar mass loss 
rate a* oc i~ 1,3 . Ciotti et al. (1991) have shown that 
galactic winds can occur during the early evolution if 
p > 1.3 and propose that the dynamically unsteady 
transition from winds to cooling flows at the current 
time may account for the strong spread observed in 
L x /Lb- However, the standard solution in Figure 2 
indicates a deficiency of interstellar gas which would 
be further reduced by galactic winds. Therefore, we 
explore here the heating effects on the interstellar gas 
using larger SNu(i„) which increases the supernova 
rate equally at all times while maintaing a cooling 
flow throughout the galactic evolution. 

The solid lines in Figure 7 show the status of the 
hot interstellar gas in our model SN1 for NGC 4472 
using SNu(t„) = 0.25, about six times larger than 
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the standard model in Figure 2. The gas temperature 
throughout the galaxy has in fact been increased, but 
only by about 20 percent, far less than the factor that 
SNu(i„) is increased. The gas density profile in Figure 
7 is also significantly improved, but is still steeper 
than the observed slope. As expected, the mean iron 
abundance in the interstellar gas (zfc) ~ 4.41 (Table 
1) far exceeds any of the values observed in NGC 4472 
and by a much larger factor than the standard model 
STD. 

Although increasing the supernova rate has a ben- 
eficial influence on both n(r) and T(r), SNu(i„) can- 
not be raised further without driving a galactic wind. 
Already when SNu(i„) = 0.3 much of the interstellar 
medium is moving outwards at time t n = 13 Gyrs and 
the gas density has declined substantially throughout 
the galaxy. The dashed contours in Figure 7 show the 
density and temperature structure at t n for the wind 
solution SN2 that develops when SNu(t„) = 0.45. 
The density has dropped to unacceptably low values 
and the temperature far exceeds observed values. 

4. MODELS WITH ADDITIONAL CIRCUM- 
GALACTIC GAS: MODELS EXG1 AND 
EXG2 

Having considered a variety of both modest and 
radical perturbations of our initial standard model 
for NGC 4472, all the same problems still remain: the 
slope of the density is too steep, the mass of hot gas 
at large galactic radii is too small, the temperature is 
too low by about a factor of 2, and the shape of the 
temperature profile is wrong. 

A clue to the resolution of this problem may lie in 
the total amount of hot gas observed in NGC 4472 
(Trinchieri, Fabbiano, & Canizares 1986; Irwin & 
Sarazin 1996; Brighenti & Mathews 1997a). At r = 
140 kpc, the largest galactic radius at which reliable 
gas density measurements are available for NGC 4472, 
the total mass in hot gas is about M gas = 3.5 x 10 10 
Mq , about 7 percent of the total stellar mass of the 
galaxy. Using our standard mass loss rate a*(t) we 
expect that approximately 8.5 percent of the total 
stellar mass of the galaxy should have been ejected 
into the interstellar medium from 1 to 13 Gyrs. How- 
ever, allowance must be made for the large fraction of 
this gas that cools. The mass budget of our standard 
model in Table 1 indicates that about 92 percent of 
the gas ejected by stars over this cosmic time interval 
has already cooled and is no longer available to emit 



x-rays. Therefore, by analogy with the standard solu- 
tion, the amount of gas observed in NGC 4472 out to 
140 kpc implies a minimum initial total gas mass of 
at least 6.3 x 10 11 Mq - this is about 10 times larger 
than the mass that can be ejected by galactic stars! 
We conclude that most of the mass that contributes to 
cooling flows in galaxies like NGC 44 7® has not come 
from mass lost by galactic stars since time 1 Gyr. 

When extra gas is introduced into the models many 
new and uncertain parameters must be determined. 
Is most of the extra gas present at early times or is 
it slowly added during the Hubble time? What is 
the origin of this extra gas: is it primordial or has it 
been enriched and ejected by galactic winds from the 
massive elliptical or other nearby galaxies? What is 
the temperature and density (entropy) of this extra 
gas? Since the additional gas must be about twice as 
hot as that in the standard solution, what total mass 
of dark matter is required to bind the additional gas? 
Although we have explored some of these questions, 
we shall not attempt to answer them in detail here. 
Our more limited objective here is to present results of 
some approximate and preliminary calculations that 
nevertheless conclusively illustrate the many benefits 
of additional hot gas in understanding the nature of 
cooling flows in bright elliptical galaxies. 

Our provisional astrophysical interpretation of the 
source of the extra gas required, as discussed below, 
is that all bright ellipticals with T gas ~ 1.5 — 2T* 
(Davis & White 1996) were originally formed in small 
galaxy groups of the sort discussed by Mulchaey et 
al. (1996). These authors have shown that a sig- 
nificant fraction of small groups of galaxies contain 
dominant bright ellipticals and that most or all of the 
observed diffuse, x-ray emitting gas in these groups 
is usually closely associated with the central massive 
elliptical. This implies that the source of the hot gas 
is closely related to the creation of bright elliptical 
galaxies, probably by mergers and tidal mass trans- 
fers. Dynamical studies of the evolution of galaxies 
in elliptical-forming groups (Merritt 1985; Bode et 
al. 1994; Garciagomez et al. 1996; Athanassoula et 
al. 1997; Dubinski 1997) require that ellipticals were 
formed early in the dynamical evolution of groups 
when mergers and tidal interactions were most likely; 
if so, both the elliptical galaxy and the hot interstellar 
gas surrounding it are very old. With this in mind, we 
consider models in which the extra hot gas is already 
present at the time t = 1 Gyr when our model calcula- 
tions begin. We shall not discuss here the alternative 
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possibility that gas gradually flows into large elliptical 
galaxies over the Hubble time. 

To construct a variant of our standard model hav- 
ing large amounts of hot gas at t = 1 Gyr, we consider 
a (somewhat arbitrary) gas density distribution at 
this early time described by n(r) — n [l + (r/r c ) 2 ]~ 3 / 4 
where n a = 0.5 cm~ 3 and r c = 0.20 kpc, out to radius 
r = 62 kpc. At r > 62 kpc we assume n(r) = 10~ 4 
cm~ 3 out to r = 470 kpc where the gas is abruptly 
terminated. The initial gas temperature T = 1.2 x 10 7 
K and iron abundance zp e — 0.1 arc uniform in this 
extended gaseous halo. We consider two models be- 
ginning with this extra gas: in the first (EXG1) a* 
has the standard value and time dependence; in the 
second (EXG2) a* is reduced by half at all times. All 
other parameters are those of our standard solution. 

As in the previous calculations, we assume that the 
dark halo of our model for NGC 4472 extends out to 
470 kpc where stars and dark halo have a combined 
total mass of 27 x 10 12 M . This is large, but com- 
parable with the masses of small groups of galaxies 
found by Mulchaey et al. (1996). The total initial 
mass of hot gas is M hot = 1.50 x 10 12 M Q within 470 
kpc, but only M hot = 0.019 x 10 12 M Q of this gas 
is within 100 kpc which includes the stellar part of 
the galaxy. Since the total stellar mass is 0.72 x 10 12 
Mq , the initial ratio of total baryonic to dark matter 
is 0.08, acceptably close to the expected cosmologi- 
cal value. It is clear from simple dynamical consid- 
erations that the dark matter halo must extend out 
considerably beyond ROSAT observations at 140 kpc 
since if the gas and dark matter were truncated at this 
radius a rarefaction wave would move into the galaxy 
in < 1 Gyr and sharply reduce gas densities and tem- 
peratures at r < 140 kpc. Therefore our assumptions 
about the mass of initial gas and dark matter are con- 
sistent with current observations of NGC 4472. 

In Figure 8 we show temperature and density pro- 
files in the hot interstellar gas at t n = 13 Gyrs for 
two models having the initial gas configuration de- 
scribed above. The first model (EXG1), shown with 
dashed lines, is computed with other parameters iden- 
tical to the standard model. The gas density is in 
excellent agreement with the observations for radii 
r > 6 kpc = 0.7r e but steepens further in. In the 
second model (EXG2) a* is reduced by a factor of 
2 so that the galactic stars contribute less to the 
density profile within r e . This model shown with 
solid contours in Figure 8 is an even better fit to 
the overall observed density profile and its luminosity, 



L x = 49 x 10 40 ergs s 1 , is very close to that observed 
for NGC 4472, L x = 45.6 x 10 40 ergs s" 1 (Eskridgc, 
Fabbiano, & Kim 1995) when scaled to d = 17 Mpc. 
For either model, the gas density near the outermost 
observed radius (~ 140 kpc) is in excellent agreement 
with observations. The temperature profile for either 
model is also in adequate agreement with observa- 
tions. This consistency with current observations of 
NGC 4472 follows from the extended region of (con- 
stant density) hot gas beyond 62 kpc that we assume 
at time t = 1 Gyr; if only this component of the ini- 
tial hot gas is retained, the agreement of the model at 
£ = 13 Gyrs with observed n(r) and T(r) is about the 
same as that shown in Figure 8. Our objective here is 
not to cosmetically adjust the parameters of the initial 
gas distribution by trial and error until the agreement 
with the present day n(r) and T(r) is perfect. Instead 
our less ambitious objectives are (i) to demonstrate 
that very hot gas can be retained by massive ellipti- 
cals throughout the Hubble time and (ii) that this gas 
is the key to understanding present-day observations. 
In a more comprehensive model, the detailed proper- 
ties of the initial (or subsequently added) gas must 
be a natural consequence of more global cosmological 
considerations. 

5. ASTRONOMICAL IMPLICATIONS OF 
CIRCUMGALACTIC GAS 

We conclude that gas ejected from stars evolving 
since the epoch of strong galactic winds, ~ 1 Gyr, 
is insufficient to account for the density and thermal 
structure observed in elliptical cooling flows. This 
conclusion is not unique to NGC 4472. Davis & 
White (1996) have shown that ROSAT temperatures 
for bright ellipticals are systematically higher than 
mean stellar temperature by ~ 1.5 — 2. Brighcnti 
& Mathews (1997a) have noticed that the normal- 
ized radial variation of gas temperature T(r/r e ) for 
six bright ellipticals having good ROSAT and ASCA 
temperatures is remarkably similar, regardless of the 
environmental context of the galaxy. Some of these 
galaxies are dominant ellipticals in small groups or 
poor clusters (NGC 5044: David et al. 1994; NGC 
507: Kim & Fabbiano 1995; NGC 1399: Jones et 
al. 1997) while others are associated with the Virgo 
cluster (NGC 4636: Trinchieri et al. 1994; NGC 
4472: Irwin & Sarazin 1996; NGC 4649: Trinchieri 
et al. 1997). Because of the apparently universal 
thermal structure among these galaxies, we conjec- 
ture that all these bright ellipticals originally formed 
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as the primary galaxy in a small group. The unifor- 
mity in the currently observed temperature profiles 
implies that the initial evolution of these ellipticals 
must have been very similar. With time some of these 
luminous ellipticals have now entered larger clusters 
where they have retained much of the gas inherited 
from the group environment. Certainly NGC 4472 is 
consistent with this interpretation since it is the cen- 
tral galaxy in a Virgo sub-cluster. Although all three 
Virgo galaxies (NGC 4636, 4472, and 4649) have the 
same positive temperature gradient at r ^ 3r e , NGC 
4649 has noticeably less hot interstellar gas than 4472 
and 4636. This may be explained by the apparent 
proximity of NGC 4649 to the center of the Virgo clus- 
ter: hot interstellar gas at large galactic radius that 
initially contributed to the production of dT/dr > 
in the central regions of NGC 4649 has since been 
ablated away by the velocity of 4649 relative to the 
intercluster medium. The pressure of the initial outer 
gaseous halo in 4649 has now been replaced by that 
of the diffuse hot gas in Virgo. Alternatively NGC 
4649 could have been originally a subordinate member 
of a small group of galaxies when its outer halo was 
tidally transferred to the group-dominant elliptical. 
These considerations are only appropriate for mas- 
sive ellipticals for which current observational data 
for p{r) and T(r) are available. The radial interstel- 
lar structure in less luminous ellipticals may differ: 
they may have less hot gas, lower overall tempera- 
tures and and steeper density gradients, resembling 
our standard model in Figure 2. Perhaps AXAF can 
answer this question in more detail. 

How is the characteristic temperature profile T{r) 
in massive galaxies produced? For further insight 
into this, we calculated the gas dynamical flow for 
model EXG1:0, a variant of model EXG1 with no 
gas contribution from stellar mass loss or supernovae. 
After 13 Gyrs the temperature of model EXG1:0 is 
in excellent agreement with observed temperatures 
of NGC 4472 within about 8 kpc from the galactic 
center (see Figure 8). At larger radii, however, the 
gas is much hotter, rising to a maximum 2 x 10 7 K 
at 35 kpc then dropping to 1.5 x 10 7 K at 80 kpc. 
At all galactic radii shown in Figure 8, the temper- 
ature of externally provided gas in our model with 
no stellar gas T(EXG1 : 0) is much higher than 
T(STD), the gas temperature resulting from stellar 
mass loss alone (Figures 2 and 8). Evidently, the rea- 
sonably successful temperature agreement of model 
EXG1 with NGC 4472 observations is due in part 



to the mixing of hot circumgalactic gas with rela- 
tively cooler gas ejected from the stars. At any ra- 
dius r ^ 100 kpc (which includes all of the stellar 
system), T(EXG1) is an approximate thermal aver- 
age of T(STD) and T(EXG1 : 0), weighted by the 
relative masses contributed at each radius by stellar 
mass loss and inflow of circumgalactic gas. By tag- 
ging the external gas and that from the stars with 
different abundances, we have determined (at time 
t = 13 Gyrs) for model EXG1 that about 60-70 per- 
cent of the gas at r < r e = 8.57 kpc comes from 
the stars while at r « 75 kpc the stellar contribution 
has dropped to 10 percent. This fractional contribu- 
tion is consistent with the relative temperatures of 
the three models shown in Figure 8 where it is seen 
that T(EXGl : 0) > T(EXG1) > T(STD) at all 
radii, although T(EXGl : 0) is closer to T(EXGl) 
at large radii and closer to T(STD) at r Ss r e . We 
recognize that this argument is only approximate be- 
cause the flow structure u(r) in these three models is 
different. Nevertheless, it is clear that the characteris- 
tic temperature profile T(r) observed in many bright 
ellipticals including NGC 4472 - rising from T — T* 
near the origin to 1.5 - 2 T* at r ~ 3r e then becom- 
ing approximately isothermal - is a signature of the 
thermal combination of two sources of hot interstellar 
gas. 

Where does the circumgalactic gas come from and 
why is it so hot? At present we cannot be entirely 
certain about the origin of this gas, but it is likely 
that its original source was "secondary infall" or the 
reversal of the Hubble flow surrounding positive den- 
sity fluctuations that grow to become massive ellip- 
ticals or small groups. Upon arrival at the young 
galaxy, the infalling gas passes through a strong shock 
that raises its temperature to approximately the lo- 
cal virial temperature T vir « (p,m p /2k)GM to t/r in 
the halo (Bertschinger 1985). It is likely that Type 
II supernovae from massive stars further heated and 
enriched the interstellar gas during the early stages of 
galactic formation. In the mass distribution for NGC 
4472 (Figure lb) T„j r increases by a factor 2.1 from 
r = r e to r = 10r e ; this may explain why the cir- 
cumgalactic gas located at r > r e is hotter than the 
stars as required in our successful models and indi- 
cated by the observations of Davis & White (1996). 
Subsequently, the hot halo gas and dark matter may 
have been dynamically transferred between galaxies 
while they were members of small groups. 

The introduction of rather large masses of hot gas 
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into ellipticals at early times may either resolve or 
complicate other problems that have received much 
attention over the years. For example, the wide scat- 
ter in L x for given optical luminosity Lb may result 
from tidal exchanges of this extended halo gas com- 
ponent between ellipticals having similar Lb within 
the small groups in which they were formed (Math- 
ews & Brighenti 1997). The general circularity of 
x-ray isophotes in rotating, massive ellipticals must 
also be understood in the context of circumgalactic 
gas. Brighenti & Mathews (1996) have shown that 
even the modest rotation of the stellar component 
observed in massive ellipticals is sufficient to produce 
strongly rotating disk configurations in the hot inter- 
stellar gas, assuming that the interstellar gas derives 
only from stellar mass loss. The x-ray images are 
flattened perpendicular to the axis of rotation. How- 
ever, if the gas expelled from stars is overwhelmed 
by inflowing low angular momentum hot gas from the 
outer galaxy, the x-ray isophotes will appear more 
circular. But it is unclear if the circumgalactic gas 
would be expected to have less rotation. [Flattened 
x-ray isophotes may nevertheless still be present in 
rapidly rotating ellipticals of moderate or low lumi- 
nosity (Brighenti & Mathews 1997b)]. The presence 
of additional gas in bright ellipticals both complicates 
and clarifies the interpretation of the iron abundance 
observed in the hot gas. Even when no circumgalac- 
tic gas is present, gas dynamical cooling flow models 
predict that the interstellar iron abundance should 
be less than that in the stars, provided the Type la 
supernova rate is very low; this occurs because the 
hot gas flows inward through a stellar system having 
a negative radial abundance gradient. But the dis- 
parity between high stellar and low gas abundances 
is so large that some have questioned the validity of 
abundance determinations from x-ray spectra (Ari- 
moto et al. 1997; Renzini 1997). The Type la su- 
pernova rate in our standard model STD is compa- 
rable with the low rate observed, but still may pro- 
duce too much iron. Nevertheless, we expect that 
observed iron abundances represent a combination of 
iron from stellar ejecta and circumstellar gas - of un- 
known abundance - that flows in from large galactic 
radii. The iron abundance in our model EXG2 is in 
fact lower than that of the standard model (Table 1), 
but not by very much for the parameters we consid- 
ered. However, Loewenstein (1996) has pointed out 
that ellipticals having the largest amount of extended, 
circumgalactic gas have on average higher iron abun- 



dances, contrary to the expected trend. 

6. CONCLUSIONS 

We have attempted without success to derive the 
known density and temperature structure in NGC 
4472 using the usual equations of gas dynamics and 
standard cooling flow assumptions. Our mass model 
for NGC 4472, based on hydrostatic equilibrium, is 
fully self-consistent and is not likely to be incor- 
rect. In our attempts to bring the results of dynami- 
cal models into agreement with observations of NGC 
4472, we considered mass dropout, spatially variable 
stellar mass loss and radiative emission rates, and 
higher supernova rates. Invariably, the total amount 
of hot gas generated by stellar mass loss that remains 
at the present time is insufficient, particularly at large 
galactic radii, and the density slope is too steep. The 
steepness of dp/dr is obviously correlated with the 
relatively low gas temperatures T w T* that result 
from thermalization of stellar mass loss and super- 
nova rates limited by the known iron abundance in 
the interstellar gas. 

All these difficulties are resolved when additional 
hot gas is provided from the outer regions of the galac- 
tic potential. We are not the first to consider addi- 
tional gas flowing into galactic cooling flows. Within 
the context of steady state solutions, both Thomas 
(1986) and Bertin & Toniazzo (1995) discussed the 
advantages of inflowing gas at the outer parts of the 
galaxy and both studies constrained their theoretical 
results to match observed properties of NGC 4472. 
The need for large gas masses in addition to that pro- 
vided by stellar mass loss is already evident from the 
high gas mass indicated by ROSAT observations of 
NGC 4472 and other similar bright ellipticals. The 
gas density in NGC 4472 at ~ 140 kpc, the outer 
limit of reliable x-ray surface brightness observations, 
is unstable to rapid outflow unless both the gas and 
the dark halo extend to significantly larger radii. We 
have therefore considered an initial configuration for 
NGC 4472 in which most of the original baryonic mat- 
ter is in an extended hot gas phase and have shown 
that much of this gas can be retained over the Hubble 
time. Our model for the initial gas configuration is 
by no means unique, but serves to demonstrate the 
possible retention of very old hot gas in massive ellip- 
ticals. 

The total mass of initial hot gas we consider within 
the large galactic halo (of radius r = 470 kpc) at 
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t = 1 Gyr, M hot = 1.50 x 10 12 M , is about 2.1 times 
larger than the stellar mass. The temperature is ini- 
tially uniform and hot, T = 1.2 x 10 7 K. This gas 
and the dark galactic halo are assumed to extend out 
to r = 470 kpc. After following the gas dynamical 
evolution, by time t n = 13 Gyr the mass of hot gas 
within r = 470 kpc is reduced by cooling and outflow 
to 0.40 x 10 12 Mq. The mass of hot and cooled gas 
within the central, stellar part of the galaxy, r < 100 
kpc, M hot = 0.04 x 10 12 M Q and M cold = 0.03 x 10 12 
M Q , are both very much less than the total stellar 
mass, M* t = 0.72 x 10 12 M Q . If the cooled gas formed 
into relatively low mass stars, the stellar mass to light 
ratio would not be affected. As gas from the reservoir 
of hot gas in the outer galaxy flows inward toward 
the stellar component, its temperature remains high, 
similar to values observed, T <~ 1.5 — 2T» and positive 
temperature gradients are generated within a few r e , 
just like the observed pattern. Although both gas and 
stars reside in the same potential, the stars are con- 
fined to just the central regions and therefore have a 
mean temperature that is less than the more extensive 
gas. Since the gas is hotter, the gas density profile is 
flattened at all radii in conformity with observation. 

The addition of hot gas to the outer regions of 
elliptical galaxies, either initially or over time, does 
not account only for the observations of NGC 4472, 
but is essential to understand similar model deficien- 
cies in all massive ellipticals. Moreover, the simi- 
larity in observed T(r) in most or all bright ellipti- 
cals observed so far (Brighcnti & Mathews 1997a) - 
both in dominant group galaxies and in members of 
rich clusters - strongly suggests that all large ellipti- 
cals have retained part of the hot circumgalactic gas 
from their previous small group environment. The 
elliptical-dominated groups observed by Mulchaey et 
al. (1996) may be typical of the early development 
of most or all giant ellipticals. Some gas dynamical 
memory of these massive gaseous halos persists even 
after bright ellipticals enter cluster environments. As 
a bright galaxy moves deep into a large cluster, the 
pressure of hot gas in its outer atmosphere may be 
supplanted by the pressure of the hot diffuse cluster 
gas. In this manner the characteristic high tempera- 
ture and positive dT j dr that are established within a 
few r e may not be greatly altered even when gas in 
the outer galaxy is stripped away; NGC 4649, deep in 
Virgo, may be an example of this. 

In the future it will be necessary to show that mas- 
sive, extended halos of hot gas, similar to the one we 



require here, are an inevitable consequence of early 
galactic evolution. The prevailing consensus today 
is that a significant outflow of gas occurred at early 
times during galactic evolution, t < 1 Gyr. The high 
metallicity observed in hot diffuse cluster gas indi- 
cates that winds were indeed common in the earliest 
times, either from the large ellipticals or from the hi- 
erarchical subunits that later merged to form them. 
It is therefore unclear how to reconcile the existence 
of the large halo of hot gas we hypothesize here with 
these early galactic winds. The hot gas is likely to 
have been acquired by massive galaxies over a few 
or many Gyrs, not provided all at once as we have 
assumed. However, our initial investigations of mod- 
els in which additional gas is supplied gradually over 
the Hubble time have not been encouraging; in gen- 
eral these models are much more complicated and less 
successful than the one we discussed here. 

But none of these complications should detract 
from our main conclusion: inflowing circumgalac- 
tic gas may help resolve many difficulties that have 
plagued the subject of cooling flows in elliptical galax- 
ies. The radial variation of gas temperature and 
density can now be understood. The long-standing 
problem of the enormous scatter in the L x /Lb plot 
(e.g. Eskridge et al. 1995) may result from dispro- 
portionate allocations of hot halo gas in group ellip- 
ticals of similar Lb following tidal transfers of halo 
material from subordinate to group-dominant ellip- 
ticals (Mathews & Brighcnti 1997). We may also 
understand why the ratio of iron abundance in the 
gas to that in the stars is so variable and often so 
small. Differing amounts of incoming circumgalactic 
hot gas, assumed to have low iron abundance, can 
dilute the iron produced by galactic stars and super- 
novae in varying amounts, leading to a variety of (of- 
ten very low) hot gas iron abundances for galaxies 
of similar Lb- Unfortunately, this simple idea is not 
supported by the few observations that currently ex- 
ist. Finally, the apparent circularity of x-ray images 
of large, rotating ellipticals may be difficult to under- 
stand without additional sources of gas since gas lost 
from the stars naturally produces easily visible x-ray 
disks out to ~ r e even in slowly rotating ellipticals 
(Brighcnti & Mathews 1996). If circumgalactic gas 
flowing in from the outer parts of massive ellipticals 
has less specific angular momentum than the stars, 
the visibility of x-ray disks would be reduced. 
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Fig. 1. — (a) Observations and our analytic fits to gas 
temperature (top) and density (bottom) for Einstein 
HRI (filled circles) and ROSAT (open circles) data; 
(b) solid curve: M tot ; long-dashed curves: stellar mass 
M* of NGC 4472; dot-dashed curves: total mass of 
hot gas; (c) solid curve: total mass density p to u long- 
dashed curves: stellar mass density p* with a hori- 
zontal cut at the core or break radius r& = 0.023r e . 

Fig. 2. — top: Filled circles and open circles are gas 
densities observed in NGC 4472 with Einstein HRI 
and ROSAT respectively. The line is the density 
variation determined by the standard model (Model 
STD) at time t = 13 Gyr. bottom: Temperatures ob- 
served by ROSAT are shown with errorbars and the 
radial dependence of gas temperature determined for 
the standard model at time t — 13 Gyr is shown as a 
solid line. 

Fig. 3. — Plots of the dropout functions A (T) (solid 
curve) and Ai(T) (dashed curve) for the ROSAT 
bandpass. 



NUA) top: gas density; bottom: gas temperature. 
The density and temperature for NGC 4472 are shown 
for comparison. 

Fig. 6. — Variation of density and temperature in 
the cooling flow model with non-uniform A (Model 
NUL) top: gas density; bottom: gas temperature. The 
density and temperature for NGC 4472 are shown for 
comparison. 

Fig. 7. — Variation of density and temperature in 
cooling flow models with increased supernova rates 
(Models SN1 and SN2). top: gas density for model 
SN1 (SNu(t n ) = 0.25) (solid line) and for model SN2 
(SNu(t„) = 0.45) (dashed line); bottom: gas tempera- 
ture for model SN1 (SNu(i„) = 0.25) (solid line) and 
for model SN2 (SNu(t n ) = 0.45) (dashed line). The 
density and temperature for NGC 4472 are shown for 
comparison. 

Fig. 8. — Variation of density and temperature in 
cooling flow models with additional gas at early times 
(Models EXG1 and EXG2). top: gas density for 
model EXG1 (normal a*) (dashed line), for model 
EXG2 (a* reduced by 2) (solid line), for model STD 
(dotted line) and for model EXG1:0 with no stellar 
gas contribution (dot-dashed line); bottom: gas tem- 
perature for model EXG1 (normal a*) (dashed line), 
for model EXG2 (a* reduced by 2) (solid line), for 
model EXG2 (a* reduced by 2) (solid line), for model 
STD (dotted line) and for model EXG1:0 (dot-dashed 
line) . The density and temperature for NGC 4472 are 
shown for comparison. 



Fig. 4. — Variation of density and temperature in the 
cooling flow model with mass dropout (Model DO) 
top: gas density; bottom: background gas tempera- 
ture of remaining gas (light solid line) and the mean 
effective temperature including that of gas dropping 
out (heavy solid line). The density and temperature 
for NGC 4472 are shown for comparison. 

Fig. 5. — Variation of density and temperature in 
the cooling flow model with non-uniform a* (Model 



This 2-column preprint was prepared with the AAS IATfrjX 
macros v4.0. 
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Table 1 

GLOBAL PROPERTIES OF COOLING FLOW MODELS 



Model 


Lx 3, 


<Z Fc /z Pc0 } b 


M hot a 


M cold a 




(10 4() ergs s" 1 ) 




(10 i0 M0) 


(10 10 M0) 


STD 


19.54 


1.03 


0.62 


6.58 


DO 


12.98 


1.15 


0.53 


6.7 C 


NUA 


29.13 


0.78 


0.76 


7.77 


NUL 


8.70 


0.97 


0.73 


6.38 


SNl d 


43.10 


4.41 


1.64 


5.01 


SN2 e 


7.66 x 10~ 4 


9.64 


1.32 x 10~ 2 




EXG1 


90.22 


0.58 f 


4.62 


7.34 


EXG2 g 


49.14 


0.85* 


4.25 


3.52 


NGC 4472 


45.6 h 


0.33' 


3.5 





a Within 100 kpc. 

b Averaged over x-ray emissivity. 

c Cold gas is distributed throughout inner galaxy. 

d SNu(t„) = 0.25 SNu. 

c SNu(t„) = 0.45 SNu. 

f Initial gas has zp e = 0.1zp e Q. 

SRor EXG2 a* is reduced by 2. 

h From Eskridge et al. (1995) scaled to D = 17 Mpc. 
'From Arimoto ct al. (1997); see text for other values. 
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